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Abstract 

We compute three- loop matching corrections to the Wilson coefficients CV and Cg 
in the Two Higgs Doublet Model by applying expansions for small, intermediate and 
large charged Higgs boson masses. The results are used to evaluate the branching 
ratio of B — > X s -y to next-to-next-to leading order accuracy, and to determine an 
updated lower limit on the charged Higgs boson mass. We find M H + > 380 GeV at 
95% confidence level when the recently completed BABAR data analysis is taken 
into account. Our results for the charged Higgs contribution to the branching ratio 
exhibit considerably weaker sensitivity to the matching scale [1q, as compared to 
previous calculations. 



1 Introduction 



In view of missing (to date) New Physics signals at the Large Hadron Collider (LHC), it is 
of utmost importance to exploit precision calculations together with precise experimental 
results in order to look for deviations from the Standard Model (SM). In this context, 
the rare decay B — > Xg'-f constitutes one of the most important processes. It is a loop- 
generated Flavour-Changing-Neutral-Current (FCNC) transition, which makes it very 
sensitive to contributions from beyond-SM particles. Moreover, its branching ratio can 
be measured with an uncertainty of a few percent and, at the same time, the result can 
be predicted within perturbation theory with a similar uncertainty. 



The current average of the measurements by CLEO [I], BELLE [2|3] and BABAR [4"f5f6] 
reads [7] 

B{B -> X s7 )\ Ei>1 . 6 GcV = (3.37 ± 0.23) x 10" 4 . (1) 

It includes the recently updated BABAR data analysis [HE]. The measurements have 
been performed with various photon energy cutoffs Eq ranging from 1.7 to 2.0 GeV. Their 
average in Eq. ([1]) involves an extrapolation to Eq = 1.6 GeV. It can be confronted with the 
SM prediction based on the Next-to-Next-to-Leading Order (NNLO) QCD calculations 
which reads B(B -> X a 7)| Sy >i. 6G eV = (3-15 ± 0.23) x 10" 4 [8p] . 

In this paper, we consider extensions of the SM Higgs sector by a second Higgs doublet, 
namely the so-called Two Higgs Doublet Models (2HDMs). They are constructed in 
such a way that no FCNC occur at the tree level (TO]- Such models have five physical 
scalar degrees of freedom, among which there is a charged Higgs boson H^ 1 that plays 
an important role for B — > X s j. We shall consider two versions of the model, usually 
denoted by 2HDM type-I and type-II where, respectively, either the same or two different 
Higgs doublet fields couple to the up- and down-type quarks. In these models, both Higgs 
doublets acquire vacuum expectation values v^ 2 such that v = \Jv \ + v \ ~ 246 GeV 
determines the W ± and Z boson masses in the same way as in the SM. The ratio V2/V1 
is denoted by tan j3. 

Comparison of the experimental results for B(B — > X s '~f) to predictions within the 2HDM 
type-II leads to the strongest constraint on the charged Higgs boson mass for tan /3 € [1, 25] 
(see, e.g., Sec. 5 of Ref. [TTJ ; the precise range depends on the treatment of uncertainties). 
So far, the constraint has been derived using only the Next-to-Leading Order (NLO) ex- 
pressions for the 2HDM contributions at the electroweak scale, while the SM contributions 
were treated at the NNLO [SJ. In the present paper, we compute the missing two- and 
three-loop NNLO matching coefficients in the 2HDM, and re-do the analysis to extract a 
lower bound on M#+ from the new experimental average (pQ). 

The outline of the paper is as follows: In the next Section, we discuss the matching 
coefficients up to the three-loop order. Besides considering the 2HDM, we also re-compute 
the SM matching contribution, and improve the three-loop results for Cj and C$- In 
Section [31 we use the new results to evaluate B(B — > X s j) to the NNLO accuracy, and 
to extract a lower bound on M H + . Section H] contains our conclusions. 

2 Matching coefficients 
2.1 Formalism 

The formalism to compute the Wilson coefficients in the 2HDM can be taken over from 
the SM case [i~2"fl~3] . We shall follow the regularization and renormalization conventions 
of those papers. In particular, we adopt the effective Lagrangian and the definition of the 



2 



operators Pj (j = 1, . . . , 8, 11) from Eqs. (2.1) and (2.2) of Ref. [13]. One should note 
that the dipole operators P-? and P§ are normalized there with inverse powers of the QCD 
coupling constant. 

Since the additional degrees of freedom of the 2HDM are all heavy, they only influence 
the Wilson coefficients Cj of the operators Pi. In order to incorporate the 2HDM contri- 
bution in a manner that is analogous to the SM analysis of Ref. [13] , we split the Wilson 
coefficients as 

qQ _ gQ,SM + £,Q,2HDM ^) 

where Q = c,t marks contributions from loop diagrams with virtual charm and top quarks, 
respectively. 

Our matching calculation is performed at the renormalization scale fiQ. It is chosen to 
be of the same order of magnitude as masses of the particles that are being decoupled 
(m t , M\y and M H +). For the NLO calculations within the SM we refer to Ref. [H]. The 
NNLO SM contributions to Cf^' SM (/xo) have been computed in Refs. p^2]fT3] . Suppression 
by m 2 c /My/ makes c^ mDM negligible. In the following, we shall consider cf mDM only. It 
is convenient to decompose them as follows 

Cf mM = ^ (0) + ^f (1) + (^) 2 ^f (2) +(^) 3 ^f (3) + .-., (3) 

where Cf is obtained from n-loop diagrams. For i = 7, 8, the tree- level coefficients 
Cf^ vanish, while Cf^ l \^o) and Cf^ 2 \fi ) have been found in Refs. [T5]I16|I17|[T8 J . In 
the present paper, we compute the three- loop corrections C^ 3 \fio) and Cg^ 3 \no), which 
constitutes the last missing element of the 2HDM Wilson coefficient evaluation that mat- 
ters for B(B — > X s ^f) at the NNLO. Furthermore, we reproduce the two-loop corrections 
Cf^ 2 \iJ, ) for i = 3,4,5,6 that have been originally found in Ref. [IH] and belong to the 
necessary NNLO matching in the 2HDM. 

In the following two subsections, we discuss our results for the Wilson coefficients in the 
SM and 2HDM. Several variables turn out to be of convenience in this context: 

M W - 2 m2 tM - , 1 1 (A\ 

w = 1 - y , r = —— — , u = 1 - r , u = 1 - - . (4) 



2.2 Standard Model 



Let us begin with recalling three-loop corrections to the matching coefficients in the SM. 
They depend on two masses: My/ and m t . In Ref. [13], expansions of the three-loop 
diagrams for My/ -C m t and My/ w m t have been performed. The final numerical results 
of that paper 

(fi = m t ) = 12.05 ±0.05, 
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y = M w I m , y = M w / m , 

Figure 1: Three- loop SM coefficients Cy (//o) and C^(jj,q) as functions of y — M w jmt 
for fi = m t . Dotted lines show their expansions around y = up to 0(y 8 ). Solid lines 
show the expansions around y — 1 and include corrections up to w 16 = (1 — y 2 ) 16 . Dashed 
lines show lower orders in w. The (yellow) band represents the physically allowed region 
for y. 



Cf\fi = m t ) = -1.2 ±0.1, (5) 

have been obtained for y = 0.488 ± 0.015. They relied on convergence behaviour and 
agreement of the two expansions for intermediate values of Mw and m t . 

wo\ t(3} 

The dotted curves in Fig. [1] show expansions of C 7 (/io = mt) and C 8 (fj,o = mt) around 
y = up to 0(y s ). Results of the expansion around y — 1 are shown as dashed and solid 
lines. The thick dashed lines include terms up to w 8 that have been provided in Ref. [T3] . 
In the present calculation, we have added eight more expansion terms for Mw m m t . 
They are shown in Fig. [1] as thin dashed lines and a solid curve that includes terms up to 
w 16 . 

t(3) 

In the case of C 7 , we observe an overlap of the two expansions for 0.2 < y < 0.35, 
which gives us confidence that the exact curve is approximated with high accuracy by the 
Taylor expansion around w = on one side, and by the asymptotic large-m 4 expansion on 
the other. The situation for C 8 in Fig. QJb) is only slightly worse. We still observe an 
improvement w.r.t. Ref. [13] due to the additional terms in the w expansion. The vertical 
bands in Fig. [1] correspond to the current experimentally allowed region y = 0.492 ±0.003 
for fiQ = m t . In the range 0.4 < y < 0.6, our improved results are very well approximated 
(to better than 0.1%) by 

C? 3) (/io = rn t ) = H-92 + 0.751 y- 1.03 y 2 , 

Cg (3) (/i = mt) = -0.764 - 2.06 y + 2.35 y 2 , (6) 

which is consistent with Eq. (JSJ) but much more accurate, and allows to substitute the 
updated value of y. In effect, the uncertainties get reduced by almost an order of magni- 
tude. 
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As far as contributions from loops involving the charm quark are concerned, the corre- 
sponding coefficients in the range 0.4 < y < 0.6 can be obtained with high precision from 
the expressions given already in Ref. [13], namely 

Cf\^ = M w ) = 1.458 2T ' 0338 , 

Cf\^ = M w ) = 1.71M/ °' 0598 . (7) 

They already provide a high-accuracy approximation in the physical region, so there is 
no need to consider higher-order terms in the expansions. 



2.3 Two Higgs Doublet Model 

In order to specify the notation, we provide the Lagrange density which defines interac- 
tions of the charged Higgs boson with fermions. Adopting the conventions from Ref. [15] , 
we have 

3 

C H+ = (2V2G F ) 1 / 2 J2 K d A u^uy i3 PL-A d m d V tJ P R )d :j H + + h.c. J (8) 

where Pl/r = (1 T 7s)/2, are the Cabibbo-Kobayashi-Maskawa matrix elements, it, 
and dj are the up- and down-type quarks with masses m Ui and m^., and Gp is the Fermi 
constant. For the 2HDM of type-I and II, the couplings A d and A u take the values 



A u = A d = —— (9) 
tan p 



and 



A u = -^- = -^—, (10) 
A d tan p 

respectively. 

The 2HDM contributions to the Wilson coefficients in Eq. ([2]) are proportional to AiA*-. 
Since the terms involving A d are suppressed by the strange-quark mass, we can safely 
neglect them. The remaining terms can be split into two parts as follows: 

-~(i,2HDM a a± ^t^HDM , a a* /it,2HDM /-, 1 \ 

For the computation of (j^ mBM ^ we can use gq_ (5.1) of Ref. [33] that has been derived 
in the context of the SM. Its application to the 2HDM is straightforward after taking into 
account that, apart from the Wilson coefficients, also the electroweak counterterms (cf. 
Eqs. (4.7) to (4.10) of p3]) and the quantities B 7 and B 8 (cf. Eqs. (3.22) and (3.23) of [T3]) 
receive additional contributions originating from the charged Higgs boson exchange. De- 
composing each of these quantities as X = X su + x 2HDM where X denotes the result 



5 




(d) (e) (f) 




Figure 2: Sample Feynman diagrams contributing to CV [(a)-(c)] and C$ [(d)-(f)] at one-, 
two- and three-loop order. 



given in Ref. [13], we obtain their 2HDM parts in D = 4 — 2e dimensions by a simple 
one-loop calculation 



7 t,2HDM 
J 2,sb 



7 t,2HDM 
J 0,sb 



AuAl T{e) 



1 2r - 1 , f . 3r - 1 m , 
+ ~ — {r e - 1) - e — — - + O (e 2 ) 



2 r- 1 



4(r - 1) 



^-2 A d A u ^ ,sb |M w ^Af H+ , 



D 



2HDM 



* D SMi 



\M W -¥M H+ i 



for 



(12) 



Both in the SM and 2HDM, the renormalization constants Z® sb and Z® sb enter the elec- 
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troweak counterterm Lagrangian in the same waj0 

CLtcr = [v; s V cb (jjP) SPr {iZ& - Zl sb m b ) b 

+ V*V tb f^^j SP R (iZ{ sb JP - Zl sb m b ) b 



(13) 



All the remaining quantities appearing in Eq. (5.1) of Ref. [13] are precisely the same as 
in the SM. 

In analogy to the SM, we have to consider vacuum integrals with two mass scales {M H + 
and m t ) in our matching calculation. Sample diagrams up to three-loops are shown in 
Fig. HJ At the one- and two-loop levels, the calculation can be performed exactly, and 
one obtains C 7 and C 8 as functions of m t /M H + [T5"fl6|ll7|ng] . At the three- loop level, 
we proceed as in Ref. [13], considering expansions around m t « Mh+, for m t <C Mjj+, 
and for m t ^> M#+. In the first case, a simple Taylor expansion is sufficient, and we 
have computed the first 16 terms in u. Calculations involving strong hierarchies require 
non-trivial asymptotic expansions. In these cases, five terms in r and 1/r have been 
evaluated. 

For the purpose of the present analysis, we have re-evaluated the Leading Order (LO) 
and NLO contributions to the renormalized Wilson coefficients, confirming the results of 
Refs. [T5|16fl7|fT8] and extending them to include higher powers in e. Such an extension 
has been necessary for renormalization at the three-loop level. Two-loop (NNLO) Wilson 
coefficients of the four-quark-operators were calculated in Ref. [19] for the MSSM. We 
have performed an independent calculation of the charged Higgs boson contribution. Full 
agreement has been found. We refrain from displaying here explicit analytical results for 
the one- and two-loop Wilson coefficients, and refer to the Mathematica file available from 
Ref. [20J. Let us only note that we set m t = mt(no) in all those lower-order terms, which 
specifies our conventions for the NNLO expressions below. 

At the three-loop level, we have been able to recover the dependence of C 7 and C 8 on /i 
by applying renormalization group techniques to the analytical one- and two-loop results. 
We find 



^7,a u a*{^) - ^?,A u At (^o — nh) 



+ ln 



00 



r (67930r 4 - 470095r 3 + 1358478r 2 - 700243r + 54970) 



in 



2187(r - 1) 
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r (10422r 4 - 84390r 3 + 322801r 2 - 146588r + 1435) , 

; 77 lnr 

729(r - l) 6 

2r 2 (260r 3 - 1515r 2 + 3757r - 1446) / 1 



Lio 1 



27(r-l) 5 V r 



1 In the corresponding Eq. (4.6) of Ref. [13], a factor of —^-e 1(i is missing in the global normalization, 
which we correct here. For clarity, the chirality projectors Pr are now displayed explicitly 
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a 



H(3) 
7AdA, 



(>o) 



a 



H(3) ^ 
,A U A* 



Oo) 



+ ln 2 



f4 

2 

m; 



-518r 4 + 3665r 3 - 17397r 2 + 3767r + 1843) 
162(r - l) 5 



r 2 (-63r 3 + 532r 2 + 2089r - 1118) 



27(r - l) 6 



lnr 



(14) 



Hnl^ 



m 



(3790r 3 - 22511r 2 + 53614r - 21069) 
81(r - l) 4 



2r (-1266r 3 + 7642r 2 - 21467r + 8179) , 
H — — -lnr 



81 (r - l) 5 
8r (139r 3 - 612r 2 + 1103r - 342) . 



Li 2 1 - - 

r 



+ ln 2 



2/^0 



m; 



27(r - l) 4 
(284r 3 - 1435r 2 + 4304r - 1425) 
27(r- l) 4 



2r (63r 3 - 397r 2 - 970r + 440) 



27(r - l) 5 



lnr 



qH{3) ^ 
8, A u A* 

+ ln 



(fjbo = m t ) 



(51948r 4 - 233781r 3 + 48634r 2 - 698693r + 2452) 



1944(r - l) 6 

r (522347r 4 - 2423255r 3 + 2706021r 2 - 5930609r + 148856) 



(15) 



lnr 



11664(r - If 
(481r 3 - 1950r 2 + 1523r - 2550) 



18(r - l) 5 



Li 2 1 



+ ln 2 (^| 

m; 



r (-259r 4 + 1117r 3 + 2925r 2 + 28411r + 2366) 



216(r - I) 1 



r 2 (139r 2 + 2938r + 2683) 



36(r - l) e 



lnr 



(16) 



a 



H{3) 
S,A d A: 



(fi = mt) + In —2 



/'n 



m 



(1463r 3 - 5794r 2 + 5543r - 15036) 



27(r - 1)^ 



r (-1887r 3 + 7115r 2 + 2519r + 19901) , 
H — — — -lnr 



54(r - l) 5 
r (-629r 3 + 2178r 2 - 1729r + 2196) . 



18(r - 1)' 



Li 2 1 - - 

r 



+ ln 2 



/i 2 , \ \r (259r 3 - 947r 2 - 251r - 5973) 



in; 



36(r- 1)' 



r (139r 2 + 2134r + 1183) 



lnr 



(17) 



18(r - l) 5 

Our results for Ci'(hq = m t ) and C^ 3 \/j,o = m t ) in terms of expansions are quite 
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lengthy. Exact values of the expansion coefficients can be found in Ref. [20]. Here, we 
present their approximate numerical values only. Considering consecutively the regions 
r — > 0, r — > 1 and r — > oo, we obtain (for Uo = m t ): 

C 7,a1'k° = 0.9225rln 2 r + 4.317rlnr - 8.278r 

-20.73 r 2 In 3 r - 1 12.4 r 2 In 2 r - 396. 1 r 2 In r - 480.9 r 2 

-34.50 r 3 In 3 r - 348.2 r 3 In 2 r - 1292 r 3 In r - 1158 r 3 

-23.26 r 4 In 3 r - 541.4 r 4 In 2 r - 2540 r 4 In r - 1492 r 4 

+42.30 r 5 In 3 r - 412.4 r 5 In 2 r - 3362 r 5 In r - 823.0 r 5 + O (r 6 ) , (18) 

C 7,a!a^ = L283 - 0.7158 u - 0.3039 u 2 - 0.1549u 3 - 0.08625 u 4 - 0.05020 u 5 
-0.02970 m 6 - 0.01740 m 7 - 0.009752 m 8 - 0.004877m 9 
-0.001721 M 10 + 0.0003378 m 11 + 0.001679 m 12 + 0.002542 M 13 
+0.003083 u 14 + 0.003404 n 15 + 0.003574 u 16 + O (u 17 ) , (19) 

C*fV 1+ = 1.283 + 0.7158 m + 0.4119 m 2 + 0.2629 m 3 + 0.1825 m 4 + 0.1347 m 5 
+0. 1040 m 6 + 0.08306 u 1 + 0.06804 m 8 + 0.05688 m 9 + 0.04833 m 10 
+0.04163 u 11 + 0.03625 u 12 + 0.03188 n 13 + 0.02827 u u + 0.02525 n 15 
+0.02269 u 16 + 0(u 17 ) , (20) 

CffV°° = 3.970 - 8.753— + 15.35- - 38.12^ + 47.09^- -103.8^ 

1 lnr 1 lnr 1 / 1 \ . , 

+79.15— - 168.3— + 24.41— - 72.13— - 274.2— + O — , 21 

iyO iyQ fjnQ lyO lyO \ rjrt\) J 

C 7,a1aY° = -20.94 r In 3 r - 123.5 r In 2 r -453.5 r lnr -572.2 r 

-8.889r 2 ln 3 r - 195.7r 2 ln 2 r - 870.3r 2 lnr - 524.1r 2 

+ 19.73 r 3 In 3 r - 46.61 r 3 In 2 r - 826.2 r 3 In r + 166.7r 3 

+36.08 r 4 In 3 r + 323.2 r 4 In 2 r + 169.9 r 4 In r + 1480 r 4 

-66.63 r 5 In 3 r + 469.4 r 5 In 2 r + 1986 r 5 In r + 2828 r 5 + O (r 6 ) , (22) 

C"f d '^ l ~ = 12.82 + 1.663 m + 0.7780 m 2 + 0.3755 m 3 + 0.1581m 4 
+0.03021 u 5 - 0.04868 u 6 - 0.09864 u 7 - 0.1306 u 8 
-0.1510 u 9 - 0.1637n 10 - 0.1712 u 11 - 0.1751 u 12 
-0.1766m 13 - 0.1763m 14 - 0.1748m 15 - 0.1724m 16 + O (m 17 ) , (23) 

c t,Saz* 1+ = 12.82 - 1.663 u- 0.8852 u 2 - 0.4827 u 3 - 0.2976 u A - 0.2021 n 5 
-0.1470 n 6 - 0.1125 u 7 - 0.08931 u 8 - 0.07291 u 9 - 0.06083 n 10 
-0.05164 m 11 - 0.04446 m 12 - 0.03873 m 13 - 0.03407m 14 - 0.03023 m 15 
-0.02702 m 16 + C(m 17 ) , (24) 
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U 7,A d A* 



a 



H(3),r->0 
8, A u A* 



a 



S,A U A* 



a 



H(3),r-+1+ 
S.A U A* 



a 



ff(3),r-K» 



H(3),r^0 



a 



■>AfjA* 



a 



H(3),r->1+ 
Ad A* 



8.088 + 9.757- 



lnr 



lnr 1 
38.43— -49.32— 



12.91- 

1 lnr 1 lnr 

-78.90— + 168.4— - 24.97— + 101.1— + 194.3 



106.2 
1 



lnr 
O 



,(25) 



0.6908 r In 2 r + 3.238 r In r + 0.7437 r 
-22.98 r 2 In 3 r - 169.1 r 2 In 2 r - 602.7r 2 lnr - 805.5 r 2 
-66.32 r 3 In 3 r - 779.6 r 3 In 2 r - 3077 r 3 In r - 3357 r 3 
-143.4 r 4 In 3 r - 2244 r 4 In 2 r - 10102 r 4 In r - 9016 r 4 

26090 r 5 lnr - 19606 r 5 



-226.7 r 5 In 3 r - 5251 r 5 In 2 r 



O (r 6 ) , (26) 



1.188 - 0.4078 u - 0.2076 u 2 - 0.1265 u 3 - 0.08570 u 4 - 0.06204 u c 



-0.04689 u 6 - 
-0.01933 m 10 



-0.009678 u 



ii 



0.03652 n 7 - 0.02907 u s 



0.01605 m 11 
- 0.008293 m 



0.01345 m 12 



0.02354 m 9 
- 0.01137m 13 



15 



0.007148 u 



16 



(27) 



1.188 + 0.4078 n + 0.2002 n 2 + 0.1190 u 3 + 0.07861 u A + 0.05531 u 5 
+0.04061 m 6 + 0.03075 m 7 + 0.02386 m 8 + 0.01888 m 9 + 0.01520 m 10 
+0.01241 m 11 + 0.01026 m 12 + 0.008575 m 13 + 0.007238 m 14 
+0.006164 u 15 + 0.005290 u ie + O (n 17 ) , 



(28 



2.278 - 5.214- 



lnr 1 lnr 1 

20.02— - 39.76— + 78.58— - 66.39— 



lnr 1 lnr 1^/1 

-91.89— + 96.35— - 300.7— + 826.2— + O — 



(29) 



-19.80 r In 3 r - 174.7rln 2 r - 658.4rlnr - 929.8 r 
-31.83 r 2 In 3 r - 612.6 r 2 In 2 r - 2770 r 2 In r - 2943 r 2 
-40.68 r 3 In 3 r - 1439 r 3 In 2 r - 7906 r 3 In r - 6481 r 3 



+54.66 r 4 In 3 r - 2777 r 4 In 2 r - 17770 r 4 lnr - 11684 r 4 

+ 1003 r 5 In 3 r - 2627 r 5 In 2 r - 29962 r 5 In r - 15962 r 5 + O (r 6 ) , (30) 

-0.6110 + 1.095 m + 0.6492 M 2 + 0.4596 m 3 + 0.3569 m 4 
+0.2910 m 5 + 0.2438 m 6 + 0.2075 m 7 + 0.1785 m 8 



-0.1546 u 9 + 0.1347n 



10 



0.1177m 



n 



0.1032 u 



12 



+0.09073 m 13 + 0.07987m 14 + 0.07040 m 15 + 0.06210 m 16 + O (m 17 ) , (31) 



-0.6110 - 1.095 u - 0.4463 u 2 - 0.2568 u 3 - 0.1698 m 4 



0.1197 u c 



10 



-0.08761 u 6 - 0.06595 u 7 - 0.05079 u 8 - 0.03987 u 9 - 0.03182 u 
-0.02577m 11 - 0.02114m 12 - 0.01754m 13 - 0.01471m 14 - 0.01244m 15 
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M H in GeV M H in GeV 



(a) (b) 

Figure 3: Three-loop coefficients C^^a* (a*o = m t) (a) and Cf A 3 K* (/io = mt) (b) as 
functions of M#+ . The dashed, solid and dash-dotted lines correspond to the expansions 
for M H + — > oo, M H + « m 4 and Mff+ — >■ 0, respectively. 



-0.01062 u 16 + O (u 17 ) , (32) 
Cfi 3) ;T°° = -3.174 + 10.89- -35.42^ + 63.744- HO.7^ + 62.264 

^) d u, iy» ,yZ iyO lyO 

In r In r / 1 \ 

-71.62— - 205.7— + 476.9— - 1003— + 0[ — ) . (33) 



At this point, a comment concerning the expansions in u — 1 — m^(/io)/Mj I+ and m = 
1 — M^ + /m((/io) is in order. For M#+ f« m 4 , one can expand either in w or in m, and 
the two expansions are easily convertible into each othero As it has been observed in 
Ref. [2T], one expects better convergence of the expansion for M#+ > m t (Mjj+ < m t ) 
if the result is expressed in terms of u (u). In the following, we shall always choose the 
better-suited representation without explicitly mentioning it. 

In Figs. [3] and HJ we demonstrate that the above expansions are sufficient to obtain the 
Wilson coefficients for any M H +. In Fig. El(a), the coefficient C^jj 3 ^* (/i = rn t ) is plotted 
as a function of Mh+. The thick-dashed, solid and dash-dotted lines show the results 
for r — y 0, r — > 1 and r — > oo, respectively, including the highest available expansion 
coefficients. Convergence of the expansions is illustrated by the thin lines that describe 
lower orders in the respective expansions. One observes that the thick dash-dotted and 
solid curves overlap for M H + w 30 — 70 GeV, which suggests that good approximations 
to the (unknown) exact results are provided by the r — > oo and r — > 1 expansions for 
M H+ < 50 GeV and M H+ > 50 GeV, respectively. Similarly, for M H + « 500 - 650 GeV, 
one observes agreement between the solid and dashed curves, which justifies the use of the 
r — >■ 1 result for M#+ < 520 GeV and the r — > result above this value. In this way, we 

2 The exact expansion coefficients from Ref. [20] should be used for the conversion. 
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Figure 4: Three-loop coefficients C 8 A ^ A + (a*o = m t) (a) and C s A * A * (/xo = fn t ) (b) as 
functions of M H +. The dashed, solid and dash-dotted lines correspond to the expansions 
for M H + — > oo, M H + « mt and — > 0, respectively. 



can define C 7 A ^ A * (/xo = "it) as a piecewise function using the expansions in the various 

limits. In Fig. &b) , the corresponding results for the coefficient C^ AA „ (/x = m t ) are 
plotted showing the same features, however, with a smaller overlap of the M#+ 3> m t and 
Mh+ ~ m t curves. For the phenomenological analysis in the next Section, we define 

[ C" { x\r^oo) for M H + < 50 GeV 
C?i 3) = I C*x\r -> 1) 50 GeV < M H+ < M 7 , x , (34) 

I tfx(r^0) M H+ >M 7:X 

with M 7iAuA * = 520 GeV and M 7tAdA * = 400 GeV. 

Analogous results for C^ 3 \fio = m t ) are shown in Fig. HI We observe the same pattern 
as for C^ 3 \/j,q = mt), which leads us to define 

\ C 8,x\ r ~> oo) for M H+ < 50 GeV 
C s,x=\ 50 GeV < M H+ < M 8iX , (35) 

with M 8AuA * u = 600 GeV and M MdA * = 520 GeV. 

The Mathematica file in Ref. [20] contains the definitions (134)) and (135|) . which allows for 
convenient numerical evaluation of the 2HDM contributions to C 7 and Cs- The updated 
SM results are included there, too. 
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Figure 5: B(B -> X s >y) in the 2HDM type-II as a function of tan (3 for M H + = 400 GeV 
(left plot), and as a function of M H + for tan/3 = 50 (right plot). Dotted, dashed and 
solid lines show central values (without uncertainties) of the LO, NLO and NNLO results, 
respectively. 



3 8{B -> X sl ) in the 2HDM to NNLO 

The framework for our numerical analysis is based on Ref. [8] where explicit results for the 
effective-theory description of B(B X s j) have been provided up to the NNLO. While 
the Wilson coefficients are known in a complete manner at this order, non-BLM NNLO 
corrections to the charm-quark-mass-dependent matrix elements (on-shell amplitudes) 
have been evaluated only in the large m c limit and extrapolated to the physical region. 

Predictions within the 2HDM to be discussed below are obtained along the same algorithm 
and using the complete NNLO matching conditions from the previous Section. However, 
two-loop purely electroweak corrections to the matching are included in the SM part only, 
as they remain unknown in the 2HDM. One should keep in mind that such electroweak 
corrections and our new NNLO QCD matching ones may be of comparable size. 

In the following, we shall discuss results for the 2HDMs of type-I and II that have been 
introduced in Eqs. (jHJ) and (TID"]) . Most of the input parameters are adopted from Ref. [8], 
except for the strong coupling constant and the top quark mass for which we use the most 
up-to-date values that are given by [22|23|24] FI 

a s (M z ) = 0.1184 ±0.0014, 

M t = (173.18 ±0.56 stat ±0.75 syst ) GeV. (36) 

The corresponding MS top quark mass equals m t (m t ) = 163.5 GeV using three- loop accu- 
racy in QCD [26"f27] and neglecting electroweak effects. As in Ref. [8], our default value of 

3 Conservatively, we use "0.0014" as the uncertainty for a s instead of "0.0007" |23I24| . This is moti- 
vated by the current tension in several precision determinations of a s (see discussion in Ref. [25]). 
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Figure 6: A closer look at the 2HDM type-II results. Left: Ai3(/i )/Ai3(/i = M H +) as a 
function of fi for tan/3 = 50 and M H + = 400 GeV at the LO (dotted), NLO (dashed) and 
NNLO (solid). Dash-dotted lines correspond to the partial NNLO result which has been 
used in Ref. [S]. Right: B(B — > X s j) as a function of M H + for tan/3 = 50. Middle lines 
show the central values, while the upper and lower ones are shifted by ±ler. Solid and 
dashed lines correspond to the NNLO 2HDM and SM predictions, respectively. Dotted 
curves represent the experimental average in Eq. ([!]). 



the photon energy cut is Eq = 1.6 GeV. Furthermore, if not stated otherwise, we choose 
fiQ = 160 GeV, fib = 2.5 GeV and /i c = 1.5 GeV for the renormalization scales, where fio 
is the matching scale, fib is the scale at which on-shell matrix elements in the effective 
theory are evaluated, and fi c is the charm quark mass renormalization scale. 

In a first step, let us discuss the branching ratio dependence on tan/3. In Fig. [5(a), we 
choose M H + = 400 GeV and show B(B ->■ X s7 ) in the 2HDM type-II for 0.5 < tan /3 < 10. 
The solid curve describes the NNLO result, while the dotted and dashed ones show the LO 
and NLO central values for comparison. One observes strong dependence for tan /3 < 2 and 
a nearly tan /3-independent result for tan /3 > 2. Actually, from tan /3 = 10 to tan/3 = 50, 
the branching ratio changes only by 0.03%. In the following, tan /3 = 50 is going to be 
our default value for the type-II model; choosing tan/3 < 2 would strengthen the lower 
limit on M#+. 

In Fig. [5(b), we show B[B — > X s ^) in the same model with tan/3 = 50 as a function of 
M H +. As expected, for large values of M^+, the 2HDM result approaches the SM one 
that overlaps with the bottom frame of the plot in the NNLO case. For Mh+ = 300 GeV 
the NNLO curve overshoots the SM prediction by about 35%, while the effect decreases 
to around 2% at M H + = 2 TeV. 

The main effect of our new three-loop terms is in reducing /i -dependence of the decay rate 
and, in consequence, stabilizing the lower bound on . This is illustrated in Fig. [6(a) 
where the charged Higgs contribution to the branching ratio AB = ^hdm — Bsm is plotted 
as a function of /iq, while normalized to its own value at /iq = Mh+ = 400 GeV. Apart 
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from the LO (dotted), NLO (dashed), NNLO (solid) curves, we also present the partial 
NNLO (dash-dotted) line that corresponds to the approach of Ref. [9]. Our calculation 
differs from the latter one precisely by including the 2HDM contributions to the NNLO 
matching. One observes a clear reduction of /io-dependence when including higher order 
corrections. Whereas the partial NNLO result for the considered ratio varies by more 
than 6.6% when fi is varied in the [80 GeV, 2M#+] range, the corresponding variation 
of our present result with full NNLO matching remains below 1.6%. The overall size of 
AB(fio) amounts to around 25% of B(B — > X S 7) SM in the considered case. 

In order to determine a lower bound on M H + we follow the approach of Ref. [9], combining 
the experimental and theoretical uncertainties in quadrature. A one-sided 95% C.L. 
(99% C.L.) bound is obtained for the lowest value of M H + for which the difference between 
experimental and theoretical central values is 1.645 (2.326) times larger than the total 
uncertainty. 

The theory uncertainty consists of four contributions that we take over from Ref. [5] . Let 
us briefly comment on each of them: 

• In Ref. [8], the non-perturbative uncertainty has been estimated to ±5%, which has 
been confirmed by the detailed investigation of Ref. [28]. We adopt this uncertainty 
for all the considered values of M#+ . 

• The charm quark mass dependence of the operator matrix elements is only partly 
known. Thus, an interpolation between the large-m c results [29] and reasonable 
assumptions for B(B — > X s ^) at m c = has to be performed. This uncertainty has 
been estimated in [8] to ±3%, which we again assume to be M H + -independent. 

• The total parametric error is obtained by combining all the partial ones in quadra- 
ture@ It amounts to around 2-^3%, however, computed from scratch for each value 
of M H+ . 

• An estimate of higher order corrections is obtained by varying the renormalization 
scales in the ranges 80 GeV < /x < max{2M H +, 320 GeV}, 1.25 GeV < fi b < 5 GeV 
and 1.224 GeV < fi c < M 6 15 = 4.68 GeV. This uncertainty is about 3 ^ 4% but 
again we compute it for each value of Mh+- 

Contrary to Refs. [8|9] , we work with asymmetric uncertainties resulting from the param- 
eter and renormalization scale variation. 

In Fig. [6](b), the NNLO result for B(B — > Xg'-f) is shown as a function of M#+, together 
with an uncertainty band that is obtained by adding all the errors in quadrature. The 
dotted (M#+-independent) curves in Fig. [6](b) correspond to the experimental result in 
Eq. ([T]). For comparison, we also show the SM prediction with the corresponding uncer- 
tainty as dashed lines. 

A correlation between the phase space factor C and m c is taken into account as described in Ap- 
pendix A of Ref. [8]. 
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Figure 7: Lower bounds on M#+ at the 95% CL. as a function of the experimentally 
determined branching ratio (abscissa) and the corresponding uncertainty (ordinate). The 
current theory uncertainty has been used in panel (a), while panel (b) presents a future 
projection with assumed reduction of the theory uncertainty by a factor of two. 



From Fig. Mjo) one can extract (using the procedure described above) the following limits 
on M H + in the 2HDM type-II: 

M H + > 380 GeV at 95% CL. , 

M H+ > 289 GeV at 99% CL. . (37) 

The above bounds replace the ones of Ref. [9] (295 and 230 GeV, respectively). Our 
95% CL. limit is very close to the one presented in Ref. |7] (385 GeV) together with the 
new experimental average (Eq. (CQ)). On the other hand, it is significantly stronger than 
the one in Ref. [I] (327 GeV) that is based on the BABAR data alone. It is interesting 
to mention that when the matching scale /j,q is varied between 80 and 400 GeV, the lower 
limits vary by around 25 GeV when our new three-loop 2HDM matching contributions 
are not included. This gets reduced to around 7 GeV only after including the three-loop 
corrections, which demonstrates the stabilizing effect of the full NNLO matching. On the 
other hand, for /i = 160 GeV fixed, the correction strengthens the limit only slightly, by 
5^6 GeV. 

The contour plots in Fig. [7] show the 95% CL. lower bounds on Mh+ as functions of 
the experimentally determined branching ratio and the corresponding uncertainty. Black 
dots correspond to the result in Eq. (JT]), while the two black lines at the bottom indicate 
the projected uncertainty to be reached by Belle II and SuperB [30]. The current theory 



16 



4 



3.6 



tan jS = 1 



2HDMtypeI 



3.4 



3 



exclusion 95% CL 




2.4 



500 1000 1500 2000 

M H in GeV 



200 



400 



600 
M H in GeV 



800 



1000 



(a) 



(b) 



Figure 8: Left: B{B — > X s ^) for 2HDM type-I with tan/3 = 1. Identification of the lines 
is the same as in Fig. [6(b). Right: 95% CL. exclusion (shaded area) for 2HDM type-I in 
the plane of tan (3 and M H + . 



uncertainty has been used in Fig. [7(a), while Fig. [7(b) presents a future projection with 
assumed reduction of the theory uncertainty by a factor of two. 

Let us now shortly discuss the 2HDM of type-I. In Fig. E(a), we show the branching ratio 
for tan/3 = 1 as a function of M^+, including the uncertainties. The meaning of the 
curves is the same as in Fig. [6(b). Contrary to the type-II model, the branching ratio gets 
suppressed with respect to the SM, and the effect becomes larger for lower values of M H +, 
which implies increased discrepancy with respect to the experimental results. Thus, it is 
possible to set a lower bound on M H + . It is shown in Fig. [E(b) where the shaded area 
represents the part of the tan/3-M#+ plane that is excluded at 95% CL. Similar results 
based on the partial NNLO predictions and the previous experimental average have been 
presented in Ref. [3T] . 

4 Conclusions 

Applying the method of expansion in mass ratios, we have evaluated three-loop matching 
conditions for the dipole operators Pj and Pg i n the 2HDM. In effect, Wilson coefficients of 
all the operators that matter for B — > Xg'y in this model at the leading order in electroweak 
interactions are now known to the NNLO accuracy in QCD. This is true not only at the 
matching scale /in but also at the low-energy scale fit, because the renormalization group 
evolution is the same as in the SM [32] . 

The main effect of including the NNLO matching is a significant reduction of /in-dependence 
of the charged Higgs contribution to the B — > X s ^ branching ratio. It is particularly trans- 
parent when considering a lower bound on M H + in the type-II model. Before taking our 
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correction into account, the bound varies by around 25 GeV when //o is varied in a rea- 
sonable range [80, 400] GeV ~ [\rn t , M H +]. Including the correction reduces the variation 
by more than a factor of 3. 

With the updated experimental average, we find that the 95% C.L. (99% C.L.) lower 
limit on Mh+ amounts to 380 (289) GeV in the 2HDM type-II. This is a universal 
(tan /3-independent) bound that can only get stronger when tan /3-dependence is taken 
into account. In practice, noticeable modifications occur for tan/3 smaller than around 2. 

In the 2HDM type-I, a 95% C.L. lower limit on Mh+ from B(B — > X s ^) can be derived 
for low tan /3 only, currently for tan /3 < 2.5 when Mh+ is above the LEP bound of around 
80 GeV [23J . In this case, considerable reduction of /i -dependence is observed, too. 

With the semi- analytical results presented in this paper, constraints on the 2HDM at the 
NNLO level can easily be updated in the future, along with developments in the measure- 
ments and in calculations of the low-energy matrix elements (that are identical in the SM 
and in the 2HDM). More precise measurements are expected in a few years from Belle- 
II [33] and Super-B [31]. On the theoretical side, the main challenges are improvements 
in analyses of non-perturbative effects [28] together with perturbative calculations of the 
NNLO on-shell amplitudes beyond the large-m c limit. In the latter case, new results 
should become available soon 135 1. 
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